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Abstract. We have obtained two deep exposures of the 
p Oph cloud core region with the ROSAT High Resolu- 
tion Imager. The improved position accuracy (l"-6") with 
respect to previous recent X-ray observations {ROSAT 
PSPC, and ASCA ) allows us to remove positional am- 
biguities for the detected sources. We also cross-correlate 
the X-ray positions with IR sources found in the ISOCAM 
survey of the same region at 6.7 and 14.3 /im, in addition 
to sources (optical and IR) known from ground-based ob- 
servations, which are young stars (T Tauri stars, with and 
without circumstellar disks, and protostars). We thus ob- 
tain the best-studied sample of X-ray emitting stars in 
a star-forming region (63 X-ray sources detected, and 55 
identified). 

We find that there is no statistically significant dif- 
ference between the X-ray luminosity functions of HRI- 
detected Class II and Class III sources, i.e., T Tauri stars 
with and without disks, confirming that the contribution 
of these disks to X-ray emission (for instance by magnetic 
reconnection between the star and the disk), or to X-ray 
absorption, must be small. 

X-ray variability of T Tauri stars can be studied by 
comparing the HRI data with the previously obtained 
PSPC data, but also using the fact that some HRI ob- 
servations were done at different epochs. The resulting 
statistics show that most of the sources are variable, and 
that their variability is consistent with a solar-like (hence 
magnetic) fiare origin. 

We use the information given both by the ISOCAM 
survey and by our HRI deep exposure to study the T Tauri 
star population of the p Oph dense cores. We confirm that 
essentially all Class II and Class III sources (embedded 
T Tauri stars) are X-ray emitters, and that a strong cor- 
relation exists between their X-ray luminosity, Lx, and 
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their stellar luminosity, L^,, with Ly^jL^, ~ 10^^. Most 
of the new ISOCAM Class II sources are not detected, 
however, which we explain by the fact that their X-ray lu- 
minosities "predicted" on the basis of this correlation are 
too faint to be detected by the HRI. 

We predict that ~40 unknown faint or embedded 
Class III sources remain to be discovered in X-rays in 
the HRI/ISOCAM overlapping area, down to a limit of 
Lx ^ 3 X lO^^ergs"^. We show that the bulk of these 
unknown Class III sources should be made of low- to very 
low-mass stars (M* < 0.1-0.6 Mq). Prospects for future 
detections with XMM-Newton and Chandra are discussed. 

Key words: (Galaxy:) open clusters and associations: in- 
dividual: p Ophiuchi cluster - stars: pre-main sequence - 
stars: formation - X-rays: stars ~ infrared: stars 



1. Introduction 

The p Ophiuchi dark cloud complex is one of the nearest 
active site of low-mass star formation (see Wilking 1992 
for a review). It is composed of two main dark clouds, 
L1688 and L1689, from which filamentary dark clouds, 
called streamers, extend to the north-east over tens of par- 
sees (e.g., Loren 1989 



de Geus et al. 1990). The main star 
formation activity is observed in the westernmost dark 
cloud, L1688, which shows a rich cluster of low mass young 
stellar objects (YSO) around two dense molecular cores, 
"core A" and "core F" in the terminology of Loren (1989) 
and Loren et al. ( 19901 ). 



The distance to the molecular complex remains some- 
what controversial (see Wilking 1992), with a usually 
adopted distance d ~ 160 pc from the Sun. From Hippar- 
cos parallaxes and Tycho B-V colors of classes V and III 
stars, Knude & Il0g (1998) have detected at d = 120 pc an 



abrupt rise of the reddening as expected from a molecular 
cloud. Based on the Hipparcos positions, proper motions, 
and parallaxes, de Zeeuw et al. ( 1999 ) gives d = 145±2pc 
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for the mean distance of the Upper Scorpius OB associa- 
tion. We adopt d — 140 pc in this article, instead of 160 pc 
used in our previous work. 

From infrared (IR) observations of s tar-fo rming re- 
gions, Lada and collab orators (e.g. Lada 1987 ; Wilking, 
Lada, & Wilking 1989 , hereafter WLY) introduced an IR 
classification and distinguished different stages of evolu- 
tion of young stellar objects (YSO). This classif icatio n was 
subsequently revisited by Andre & Montmerle (1994, here- 
after AM) to incorporate results of millimeter continuum 
studies on circumstellar dust. The IR sources are classi- 
fied in three classes, according to their spectral energy 
distributions (SEDs). This classification, initially defined 
empirically, is now well understood in terms of evolution 
of low-mass stars at their earliest stages. Submillimeter 
observations led to the discovery of cold objects, younger 
than the IR sources, and thus to the introduction of a 

2000|) . 



1992 



fourth class named "Class 0" (Andre et al. 
Class sources are very young protostars, peaking in the 
submillimeter range, at the beginning of the main accre- 
tion phase. Class I sources are evolved IR protostars, opti- 
cally invisible, in the late accretion phase. Class II sources 
are YSO surrounded by optically thick circumstellar disks. 
Class III sources are YSO with an optically thin circum- 
stellar disk or no circumstellar disk. Studies of optically 
visible YSO, T Tauri stars, led to another classification 
based on the Ha line, which separates "classical" T Tauri 
stars (CTTS) from "weak-Hne" T Tauri stars (WTTS) 
according to their equivalent width in emission, with a 
boundary at EW[IIq]^ 5 — 20 A, depending on the spec- 
tral type (Martm 1997). CTTS and WTTS are usually 
taken to be identical to Class II and Class III sources re- 
spectively, on the basis of their IR SED (see AM for a 
discussion about these two classifications). We will asso- 
ciate in this article Class II (Class III) sources with CTTS 
(WTTS). 

Several ground-based near-IR surveys (e.g. Wilking 



et al. 1989 


Greene et al. 


1994 


hereafter GWAYL; Bar- 


Sony et al. 


1997, hereafter BKLT; and references therein) 



discovered in a «1 square degree area around the densest 
regions (with survey completeness limit down to ivT ~ 14) , 
~100 low-luminosity embedded sources. More recently, 
the ISOCAM camera on-board the Infrared Space Obser- 
vatory satellite imaged a half square degree centered on 
LI 688 in the mid-IR (LW2 and LW3 filters, respectively 
centered at 6.7 /zm and 14.3 /xm - ISOCAM centr al pro - 
gramme surveys by Nordh et al.; see Abergel et al. 1996|) , 
and recognized 68 new faint young stars with infrared ex- 
cess (Bontemps et al. 200C). 

Near-IR spectroscopy has been used to determine spec- 
tral types of an increasingly large number of p Oph YSO 
(see the pioneering works of Greene & Meyer 1995|, and 
Greene & Lada 1996). Recently, Luhman & Rieke (1999) 



with a representative population from the outer region of 
the cluster {K <ll). 

The p Oph dark cloud YSO have also been extensively 
studied in X-rays. Early observations with the Einstein 
Observatory satellite showed that at the T Tauri star stage 
YSO are bright and variable X-ray emitters in the 0.2- 
4keV energy band (Montmerle et al. [19831) . When the S/N 
ratio is sufficient large, their X-ray spectra can be fitted 
by a thin thermal model, with temperatures « 1 keV and 
absorption column densities Nh ^ 10^°-10^^ cm~^. Vari- 
ability studies and modeling led to explain the X-ray emis- 
sion in terms of bremsstrahlung from a hot (Tx ~ 10^ K) 
plasma trapped in very large magnetic loops, in other 
words in terms of an enhanced solar-like flare activity (see 
reviews by Montmerle et al. 1993| ; and Feigelson & Mont- 
merle 1999, hereafter FM). Casanova et al. (1995) - here- 
after CMFA - reported deep ROSAT Position Sensitive 
Proportional Counter (PSPC) imaging of the p Oph cloud 
dense cores A and F. They detected in the 35' x 35' central 
portion of the field (the inner ring of the ROSAT detec- 
tor entrance window support structure) 55 X-ray sources 
in the 1.0-2.4 keV energy band. For three X-ray sources, 
one or several Class I sources lie within the error boxes of 
X-ray peaks, but other counterparts are possible (unclas- 
sified IR sources, T Tauri stars). X-ray emission from one 
of these Class I sources, YLW15 (=IRS43 in WLY), was 
unambiguously confirmed with a follow-up ROSAT High 
Resolution Imager (HRI) observation by Grosso et al. 
( |1997| ). The outer p ortion of the CMFA PSPC field, ana- 
lyzed by Casanova (1994), contains 36 X-ray sources. The 
optical spectroscopic classification of these X-ray sources 
and other X-ray selected stars in the p Oph dark cloud 
vicinity, based on Hq and Li I (670.8 nm) spectroscopy, 
was made by Martin et al. ( 1998] ), doubling the number of 
PMS stars spectroscopically classified in the p Ophiuchi 
area. 

Observations of harder X-ray (>4 keV) from the p Oph 
dark cloud were initially only possible with non-imaging 
instruments. Tenma and Ginga revealed unresolved emis- 
sion from the cloud core region, with a hard X-ray spec- 
trum with fcTx ~ 4 keV and Nh ^ 10^^ cm~^ (Koyama 



1987 ; Koyama et al. 1992). Wide-energy band imaging ob- 
servations became possible with ASCA in the rang e 0.5- 
lOkeV. In the p Oph dark cloud, Koyama et al. (1994) 
detected hard X-rays from T Tauri stars, with fcTx up to 
~8keV in the case of the WTTS DoAr21. There is also 
some evidence for unresolved hard X-ray emission from 
embedded young stars below the point source de tectio n 
limit. From this ASCA observation, Kamata et al. (1997), 



obtained if-band spectroscopy for ~100 sources, combin- 
ing a magnitude-limited sample in the cloud core {K < 12) 



found additional T Tauri stars and detected three X-ray 
sources associated with Class I sources, but with large X- 
ray error boxes (15"-30"). 

There is a deep connection between IR and X-ray ob- 
servations of star-forming regions. Sensitive ground-based 
near-IR surveys penetrate dark clouds (except for dense 
cores) so that their source populations are frequently dom- 
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inated by ordinary stars in the Galactic disk. Space-based 
mid-IR isolates YSO with significant circumstellar mate- 
rial and effectively eliminates the background star pop- 
ulation, but they will miss the recognition of YSO with 
less massive or absent disks. X-ray emission, in contrast, 
is elevated by 1-4 orders of magnitude in YSO of all ages, 
irrespective of a disk presence. It thus provides a unique 
tool for improving the census of young star clusters. 

In this article, we present the results from the HRI 
follow-up of the CMFA PSPC observation. The high an- 
gular resolution of these observations allows us to find 
counterparts to all X-ray sources without ambiguity. The 
comparison with the sensitive ISOCAM survey of the p 
Oph dark cloud significantly improves the existing classifi- 
cation of these counterparts and allows us to do statistical 
studies on a well defined sample. 

We first present the ROSAT HRI observations: im- 
age analysis, source detection and identification (§|]). We 
incorporate the ISOCAM survey results from Bontemps 
et al. ( ^OOC ) and we present the resulting IR classifi- 
cation for the HRI sources (§||). The next sections dis- 
cuss the X-ray luminosity of the fl^fl/ detected TTS (§|), 
and the X-ray detectability of the embedded TTS popu- 
lation (§^. Next (§^, we show that the HRI census of 
Class III sources cannot be complete, and that numerous 
unknown low-luminosity Class III sources, perhaps includ- 
ing brown dwarfs, must exist. Summary of the main results 
and conclusions are presented in §0, where prospects for 
improvements with XMM-Newton and Chandra, are also 
discussed. 

Appendix A gives details about the HRI X-ray source 
detection, and lists the X-ray detections. Optical finding 
charts, and identification list of the HRI X-ray sources 
can be found in Appendix B. Appendix C compares these 
HRI observations with previous PSPC ones. Appendix D 
discusses the status of optical/IR counterparts without IR 
classification. 



2. The ROSAT HRI X-ray observations 

2.1. The ROSAT HRI images 

We have observed the dense cores A and F of the p Oph 
dark cloud with the ROSAT HRI. The detector is sen- 
sitive to the 0. 1-2.4 keV energy range, but has no spec- 
tral resolution. The two observation fields were respec- 
tively centered on approximatively the WTTS DoAr21 
(a = 16h26'"2!4, 5 = -24°23'24" [J2000]) and on the 
Class I protostar YLW15 = IRS43 (q = 16''27'°26?4, 5 = 
-24°40'48" [J2000]). ROSAT HRI images have a diam- 
eter of ~ 40'. Fig. |l] displays the two observation fields, 
which include the dense DCO+ cores A, B, C, E, and F 




Fig. 1. Map of the p Ophiuchi dark cloud observation 
fields. 2°diameter ROSAT PSPC field of view, a nd it s 
central 35' x 35' studied by Casanova et al. (1995); 
40'diameter ROSAT HRI fields of view and YLW15 area 



(Loren et al. 1990[ ), most of the area studied by CMFA 
and ASCA, as well as the ISOCAM survey. The first ob- 
servation field, centered approximatively on p Oph A, will 
be referred to as the "core A field" : the second observation 



published in Grosso ct al. (|1997| ); ASCA G/S* field of view 
from Koyama et al. (1994) and Kamata et al. (1997); 
ISOCAM survey field (central programme observations by 
Nordh et al.; see Abergel et al. 1996). The background im- 
age of the p Ophiuchi star forming region is taken from 
the first Digitized Sky Survey, contrast was enha nced with 
Gaussian histogram specification (Canzian |l997 ). Regions 
of high visual extinction are clearl y visible. DC0+(J=2-1) 
contours from Loren et al. ( 199C| ) show for reference the 
location of dense cores A, F. 



field, centered on dense core F, will be called the "core F 
field". The core A field was observed between 1995 Au- 
gust 29 and 1995 September 12 with a total exposure of 
51.3 ks. The core F field was observed at three different 
epochs (hereafter observations #1, #2, and #3): between 
1995 March 9 and 14 (12.5 ks), between 1995 August 18 
and 20 (27.5 ks), and between 1996 September 7 and 11 
(37.2 ks). These three different epochs give a total expo- 
sure of 77.2 ks (see Table |l| for the log of ROSAT HRI 
observations details). 

2.2. The HRI image analysis 

We have analyzed separately the four data sets. Stan- 
dard source detection algorithms were used to find X- 
ray sources, and to search optical counterparts for each 
set. We then selected X-ray sources with the best posi- 
tion accuracy (usually the brightest not too far away from 
the field center) and having an unambiguous counterpart. 
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Table 1. Log of ROSAT HRI Observations. 









Begin 


End 




Elapsed 




Pointing 


Core 


Obs. 










Time 


Exposure 
[ks] 


ID 


# 


Date 


Hour (UT) 


Date 


Hour (UT) 


[ks] 


201835h 


A 


1 


1995 Aug 29 


18''31""30" 


1995 Sep 12 


22"35"'47" 


1224.3 


51.3 


201834h 


F 


1 


1995 Mar 09 


00'"51""03" 


1995 Mar 14 


05''16'"43= 


447.9 


12.5 


201834h-l 


F 


2 


1995 Aug 18 


19''29™21= 


1995 Aug 20 


14''52"'05' 


27.8 


27.5 


201834h-2 


F 


3 


1996 Sep 07 


23'"12""46= 


1996 Sep 11 


13''46'"23" 


311.6 


37.2 




F 


1+2+3 


1995 Mar 09 


00''51™03" 


1996 Sep 11 


13'M6"23= 


47 739.3 


77.2 



corrected all X-ray positions from the existing offsets, and 
used this improved astrometry to remove possible ambi- 
guities in the identification of the X-ray sources. Details 
can be found in Appendix A. 

We find 63 HRI X-ray sources. Fig. || shows the po- 
sitions of these sources, superimposed on a combined IR 
and optical image of the p Oph cloud. Coordinates, error 
boxes, likelihoods of existence, and count rates of these 
HR I X -ray sources can be found in Appendix A (Ta- 
ble O). 



2.3. HRI X-ray source identifications 

Identification of the ffi?/ X-ray sources was made by cross- 
correlations with published lists of confirmed or suspected 
cloud members (AM), IR sur veys (Greene & Young 1992; 



BKLT; and Bontemps et al. |2000| , see g3h, if -band spec- 



troscopic (Luhman & Rieke 1999|) , radio sur veys ( Andre 



Mont merle, & Feigelson |1987| Stine et al. [1988| ; Leous 
et al. |l99l| ), and with previ ously known X-ray sources 
(PSPC: CMFA, Martin et al. |l998| ; ASCA: Kamata et al. 
1997| ). For X-ray sources_ without published counterparts 



we have used 5'/MBj4I?,[| and we have also searched op- 
tical counterparts on optical red band images from the 
Digitized Sky Survey.^] 

As shown by the finding charts in Appendix B, thanks 
to the good angular resolution of the ROSAT HRI ( PSF 
FWHM ~ 5" on axis; due to the mirrors as well as the 
detector), the position accuracy (l"-6", see Col. 6 Ta- 



ble A.l) allows us to find counterparts almost without 
ambiguity.^ Identification lists for the core A and core F 
fields are given in Appendix B , Table B . 1 . We also discuss 
in Appendix B the identification of X-ray sources with a 
low statistical significance. 

Nearly 90% of the HRI X-ray sources are identified. 
We detect only - 70% of the PSPC X-ray sources (CMFA; 



On-line version at 
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http:/ /simhad. u~ strashg.fr/ Simbad 



org/ cgi— bin/ dst 



When two possible counterparts are in the error box (this 
happens o nly th ree times, see the finding charts in Appendix B, 
and Table 



B.l 



notes), we take the more luminous in the J- 
band, since the X-ray luminosity of TTS is correl ated with the 
J-band luminosity (see CMFA, and below §U, §[^21). 



Casanova 1994), but this can be explained by the differ- 
ence in sensitivity between the two instruments, and the 
intrinsic variability of the X-ray sources (see Appendix C). 

3. Nature and IR properties of the HRI sources 

Following the results of CMFA, essentially all X-ray 
sources we found in the p Ophiuchi dark cloud should be 
young stars, a number of them being still embedded in the 
cloud. Embedded YSO are mainly studied at IR and mil- 
limeter wavelengths, but they may all be potential X-ray 
emitter regardless of their IR classif ication. Indeed, using 
the results of Wilking et al. (|l989D , AM, and GWAYL, 
CMFA analyzed their results in the light of the IR obser- 
vations of the stars they observed in X-rays. We can then 
(i) use the published IR surveys to provide a list of recog- 
nized YSO members of the cluster to be compared to the 
observed X-ray properties, and (ii) use X-ray emission to 
discriminate between true cluster members and the many 
potential background stars seen in IR images. 

3.1. New Class II and Class III source census after 
ISOCAM 

Near-IR surveys of the p Ophiuchi cluster, sensitive 
enough to detect low-luminosity embedded young stars, 
have recently been published (Comeron et al. 1993 ; Strom 
et al. |1995| - hereafter SKS; BKLT). However these 



ground-based surveys encountered limitations in recogniz- 
ing the nature of all the embedded sources, and have there- 
fore not much increased the number of bona fide members 
of the p Oph cluster. The mid-IR camera aboard ISO, 
ISOCAM, produced a ma p of the p Oph main cloud, used 
by Bontemps et al. (200C) to study the young star popu- 
lation. This mid-IR study resulted in a significantly more 
complete census of the p Oph cluster population. We here 
use this new census as a basis for discussion about the 
nature of the detected X-ray sources and to estimate the 
occurrence and properties of the X-ray emission of the 
different classes of YSO. 

The mid-IR photometry at 6.7 and 14.3 /zm appears 
invaluable to characterize sources with IR excesses, i.e., 
Class I sources and C lass I I sources (e.g., Nordh et al. 
1996 ; Bontemps et al. 1998 ). In p Oph, Bontemps et al. 



have doubled the number of Class II sources known. 
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Fig. 2. X-ray sources in the ROSAT HRI fields. The composite ISOCAM map of Abergel et al. ( |1996D (LW2 fiher 
image [5-8.5 /im], plus LW3 filter image [12-18 fim]) is merged with the background optical image taken from the first 
Digitized Sky Survey. DC0+(J=2-1) contours show the location of dense cores named A, B, C, E, F by Loren et al. 
(199C). A scale of 0.5 pc is shown for d — 140 pc. The positions of the X-ray sources are marked by 50"-size squares. 
Typical X-ray error positions range between l"-6". 



These authors conclude that the sample is complete down 
to stellar luminosities, L^,, as low as 0.03 L0, thus extend- 
ing downwards the luminosity function obtained from the 
ground by about one order of magnitude. However, just as 
in the near-IR, these measurements alone cannot charac- 
terize the nature of sources without IR excess: these can 



be either Class III sources (diskless TTS) or background 
sources. 

The tentative Class III classification of a number of 
sources coming from previous X-ray observation (CMFA), 
or from this article, is now confirmed by ISOCAM, which 
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detected no IR excess. We call these sources "new" 



Class III sources (see Table |B.l| , and |5|) . 

Since we already identified one of the HRI sour ces a s 
a foreground star (ROXF37 = HD148352, see Table 
the question arises that a number of these new sources 
with Class III spectra may also be field stars, conta minat - 
ing the genuine young star sample. Guillout et al. ( 1996 ) 
has estimated the stellar content of flux-limited X-ray sur- 
veys, based on the age-dependent stellar population model 
developed by the Besangon group. For the average sen- 
sitivity of our ffi?/ observation ('-^ 6 x 10~^ctss~^, see 
Fig. the p Oph galactic latitude (~ 20°), and our to- 
tal field of view (0.5 square degrees), this model yields 
an estimate of 5 contaminating stars (P. Guillout, private 
communication). Removing IID148352, this leaves 4 pos- 
sible field star candidates among the 21 "new" Class III 
sources and "Class II-III" sources listed in Table B.l. This 



number is small enough not to affect significantly our dis- 
cussion, and in what follows we will simply neglect the 
possible contamination of our various Class III samples 
by field stars. 

The Class III source population will be studied in de- 
tail below (Sg and §|). 

3.2. Cloud extinctions and stellar luminosities from 
near-IR photometry 

YSO suffer from large amounts of extinction by dust - 
cloud dust plus circumstellar dust - (up to Ay = 60 or 
more in p Oph), which strongly affects their fluxes at all 
wavelengths of interest, but does not necessarily prevent 
their detection in soft X-rays (see, e.g., CMFA). Despite 
this effect, it is possible to estimate the total extinctions 
along the line of sight and stellar luminosities. Near-IR 
photometry data appear to provide the most reliable es- 
timate for the luminosity of the embedded young stars 
(see discussion in Bontemps et al. 2000 and references 
therein): the J-band fluxes are usually almost purely pho- 
tospheric and thus trace the stellar luminosities very well 
(e.g., GWAYL; SKS). Similarly the J - color is a sensi- 
tive tracer of interstellar extinction. We therefore use the 
extinctions and stellar luminosities derived from near-IR 
photometry for the ISOCAM sample by Bontemps et al. 
(|2000D . 

In CMFA the bolometric luminosities were approx- 
imately estimated from the J-band fluxes using the 
GWAYL conversion based on an observational correla- 
tion between J fluxes and bolometric luminosities for Tau- 
rus and Chameleon TTS. The bolometric luminosity com- 
prises in principle the total luminosity (accretion -I- stel- 
lar) of a YSO. However for young PMS stars like Class II 
sources and Class III sources, this luminosity should co- 
incide with their stellar luminosity since their accretion 
luminosity (if any) has become significantly smaller than 
the photospheric luminosity. Finally we note that the con- 
version between the J-band absolute magnitude and the 



stellar luminosity used by Bontemps et al. (2000) is nu- 
merically similar to the CMFA conversion between the 
dereddened J-band and the bolometric luminosity. 

Cols. 8-9 of Table B.l give for each source its more 



up-to-date values of interstellar extinction and stellar lu- 
minosity, determined by Bontemps et al for 140 pc. The 
reader will fi nd the details of the calculations in Bontemps 
et al. (|2000| ). 



3.3. Nature of the HRI sources from IR data 

The YSO evolutionary stage, inferred from IR spectral 
energy distributions (see §^]^), is available for most of 
the X-ray sources. The resulting census of the 55 X-ray 
sources with stellar counterparts is: one Class I protostar 
(YLW15=IRS43); 23 Class II sources, including 4 new 
ISOCAM Class II sources; 21 Class III sources, includ- 
ing 13 new Class III sources^; 8 whose classification is 
either Class II or Class III sources (see discussion in Ap- 
pendix D); one Class III e arly-ty pe star (SI, with spectral 
type B3; see Andre et al. 198S| ), and one main sequence 
foreground star (IID148352 with spectral type F2V; see 
the Hipparcos catalogue). 

In Table p3.l| , Col. 6 lists cross-identifications with the 
survey: red (blue) means /S'OC^M sources with 
(without) IR excess. Col. 7 gives the IR classification. 
Cols. 8-9 the extinctions and stellar luminosities (from 
Bontemps et al. pOOC ). 

Therefore, the present analysis has revealed no new 
X-ray emitting Class I source, apart from the Class I pro- 
tostar YLW15 which has been the subject of a specific 
study (Grosso et al. |1997D . The HRI/ ISOCAM sources 
are overwhelmingly TTS, for which improved results are 
described in the following sections. 



4. X-ray luminosities of the HRI-detected T Tauri 
stars 

4.1. Derivation from source count rates 

We a ssume a fiducial TTS X-ray spectrum (see Mont- 
merle 1996) having fcTx=l keV plasma, with cosmic abun- 
dances, and Raymond-Smith line emissivities, and with 



i nterst ellar absorption based on Morrison & McCammon 
( 1983) cr oss sections. We use the relation given by Ry- 
ter ( 1996 ) to estimate the hydrogen column density, Nh, 
from the visual extinction. Ay, determined from IR data^ 
Nh = 2.23 X lO^Mycm^^. The intrinsic (i.e., extinc- 
tion corrected) X-ray luminosities in the full ROSAT en- 



* 4 were already PSPC Class III source candidates, and are 
confirmed by ISOCAM, which detected no IR excess; 3 were 
probably detected by the PSPC, but due to its lower angular 
resolution the o ptic al or IR counterpart was uncertain (see 
and attached notes); 



Col. 2 of Table 



B.l 



are genume new 
X-ray detections, probably resulting from variability. 
^ Noted Av.iR in CMFA. 
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ergy band (0.1-2.4 keV), Lx, were calculated for classi- 
fied sources from source count rates for d ~ 140 pc using 
EXSAS, and are given in Cols. 10-13 of Table ^ The 
intrinsic X-ray luminosities span the range Lx ~ 1 x 10^^ - 
~ 4 X 10'^^ erg s~^. (For X-ray sources without extinction 
estimates we label the detection exposure with a question 
mark.) 



^.2. Luminosity Junctions of HRL detected Class II and 
Class III sources 

We here compare the extinction corrected X-ray luminosi- 
ties of the Class II and Class III sources detected by the 
HRI'm the ZSOCMM field^ in order to evaluate the contri- 
bution of the circumstellar disk to the X-ray absorption, 
or to the X-ray emission (for instance by magnetic recon- 
nection between the star and the disk; see Montmerle et al. 
2000| ). 

Fig. H shows the cumulative X-ray luminosity distri- 
bution functions for these two populations, estimated us- 
ing the ASURV statistical software package (rev. 1.2;^ La 
Valley et al. 1992| ), which takes upper limits into account. 
These distributions are mathematically identical to the 
maximum likelihood Kaplan-Meier estimator (Feigelson & 
Nelson |1985D . Their mean X-ray luminosities (in ergs ^) 
are given by < log(Lx) >= 30.3 ±0.1 for Class II sources, 
and < log(Lx) >= 30.4±0.2 for Class III sources. We used 
nonparametric two-sample tests implemented in ASURV 
— Gehan's generalized Wilcoxon test, Logrank test, Peto 
& Peto generalized Wilcoxon test, Peto & Prentice gen- 
eralized Wilcoxon test — to see whether the difference 
between the two luminosity functions is significant. These 
tests gave a high probability (46-72 %) that they are sta- 
tistically indistinguishable. This result agrees with previ- 
ous deep studies of the p Oph main cloud (CMFA ) and 
Chamaeleon I (Feigelson et al. 1993 ; Lawson et al. 1996 ) 
YSO populations. 



In contrast, Neuhauser et al. (1995) found in the 
Taurus- Auriga star-forming region that Class III sources 
are more X-ray luminous than Class II sources. Neuhauser 



et al. ( |1995D used the ROSAT All Sky Survey (RASS) to 
cover a large area of the Taurus- Auriga (~ 900 square de- 
grees), including dense cores, but also away from this star- 
forming region. This shallow survey and large area must 
be compared to our deep pointed observations, where our 
field of interest covers 0.5 square degrees. In our observa- 
tions we focus only on the dense cores, studying a younger 
population of YSO. Contamination of the T Tauri star 
sample by a more evolved, wide-spread, and older Class III 



22 Class II sources and 19 Class III sources; for the core F 
obser vatio n we take the mean X-ray luminosity in Col. 13 of 
Table 
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Fig. 3. Cumulative normalized X-ray luminosity functions 
of Class II and Class III sources detected with the HRI in 
the p Oph cloud. The dashed (dotted) histogram shows 
the Class II (Class III) source integral X-ray luminosity 
functions. The solid histogram shows the total cumulative 
X-ray luminosity functions. The straight line shows the fit 
for the total cumulative X-ray luminosity functions. The 
left scale gives the Kaplan-Meier estimator. The right scale 
gives the cumulative number of X-ray sources. 



source population, may explain the discrepancy with the 
result of Neuhauser et al. ( 1995 ).P| 

Thus, we confirm that the contribution of the disk of 
Class II sources to their X-ray emission, or to X-ray ab- 
sorption, must be small. 

Next, we combine the X-ray luminosity distributions 
of Class II and Class III sources to obtain the cumulative 
X-ray luminosity distribution function of all TTS detected 
by the HRL For log(Lx) between 29.7 (log(Lx,break)) 
and 31.3 (log(Lx,max)), the distribution is loghnear: M{> 
log(Lx)) = -21.9 X (log(Lx) - 31.3). For log(Lx) < 
log(Lx, break), the distribution shows a downwards trend, 
which is due to our lower efficiency to detect weak X-ray 
emitting TTS. We deduce from this linear relation the to- 
tal X-ray luminosity emitted by the X-ray sources with 
X-ray luminosity between Lx,max and Lx,mm (an arbi- 
trary value): Lx,tot = 21.9 x (Lx,max - ix,min)- Thus, the 
total X-ray luminosity of this group of X-ray sources is 
dominated by the brightest sources (with ix,max), and is 
very weakly dependent on the loglinear fit. For ix.min — 
ix, break, we find Lx,tot ~ 4.3 X 10^^ ergs^^. This value is 

ergs" 



close to the asymptotic value 4.4 x 10^^ ergs ^ obtained 



B.l 



Version available at http://www.astro.psu.edu/statcode. 



^ In addition one should note that the RASS includes the 
soft band (0.1-1 keV), which is more sensitive to the extinction 
than the hard band (1-2.4 keV) taken by CMFA. 
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taking Lx,min = 0, thus future detections of new X-ray 
emitting TTS with low X-ray luminosity will not greatly 
affect this result. 



4-.3. HRI source variability 

We present in Appendix C a study of the variability of the 
X-ray sources which were observed both by the HRI and 
the PSPC, and we show that some sources were in a high 
X-ray state during the HRI or the PSPC observations. 
Here, we study the variability of the Core F HRI sources. 

Core F observations comprise three time-separated ob- 
servations, which allow us to reiterate the "Christmas 
tree" luminosity study made with Einstein Observatory 
by Montmerle et al. (1983). The idea suggested by the 



similarity between Class II sources and Class III sources 
in X-rays was to assume that all the X-ray sources are 
basically one single type of X-ray object, seen in differ- 
ent states. The result was that the distribution of the flux 
variations could be approximated by a power-law with an 
index (3 — —1.4. 

We have estimated whenever possible for each Core F 
HRI source the X-ray flux variations from the observed 
high/low luminosity ratio, TZ = -Lx,high/-^x,iow, based on 
the three observations. This yields 27 values including 11 



1.0 



a = -0.75±0.03 



10 




o 



lower limits. The integral distribution for a given ratio 
is estimated using the maximum likelihood Kaplan-Meier 
estimator (see Fig. ^. We find a power-law distribution: 
A/'(> Lx,high/ix,iow) oc (Lx,high/ix,iow)", with an index 
a = —0.75 ±0.03. This implies that the differential distri- 
bution dM /dR. follows a power-law distribution of slope 
/3 = a - 1 = -1.75. 



We suggest following Montmerle et al. (1983) that this 
power-law behavior, may be explained in terms of variabil- 
ity due to stellar flares, if interpreted in terms of stochas- 
tic relaxation phenomenon (Rosner & Vaiana 1978), and 
dominating the X-ray activity of the underlying stars. 
Such a power-law behavior is seen in the solar flares in ra- 
dio, optical, soft and hard X-ray emission with the power- 
law index -/?= 1.1-3.0 (see review in Aschwanden et al. 
1998). For soft X-ray emission — /3=1.7-1.9, which is con- 



sistent with our result, and supports the analogy with the 
solar magnetic activity. 

5. X-ray detectability of the embedded T Tauri 
star population 

We use the information given both by the ISOCAM sur- 
vey and by our HRI deep exposure to study the X-ray 
detected TTS population of the p Oph dense cores. We re- 
strict the following studies to the HRI/ ISOCAM overlap- 
ping area.0 This area comprises 98 Class II sources (clas- 
sified from ground-based and ISO observations), includ- 
ing 52 new ISOCAM Class II sources,[3 and 35 Class III 
sources (characterized as YSO from X-ray or radio obser- 
vations, and classified as Class III sources from ground- 
based and ISO observations), including 21 new Class 111 
sources {HRI or PSPC X-ray sources without IR excess 
observed by ISOCAM). We will call these sources the 
"TTS sample". 

Our HRI observation detected a large number of 
sources, yet these constitute only 30% of the "TTS sam- 
ple" . In this section, we examine the reasons why the other 
TTS were not detected, and in particular whether the un- 
detected TTS form a separate population of genuinely X- 
ray weak objects. 



0.1 L , , , , ... \ L I 

1 10 
log(l..,i„j/T.,,,..,) 

Fig. 4. Distribution of X-ray luminosity variations. We 
have estimated whenever possible for each X-ray source of 
core F the X-ray variability amplitude with the ratio high- 
on low-observed luminosity during the three observations. 
The frequency distribution has been estimated using the 
maximum likelihood Kaplan-Meier estimator. We find a 
power-law distribution: N{> ^^^f^) oc (-^yS^^)", with 
a — —0.75 ± 0.03, consistent with a variability due to 
stellar flares, which supports the analogy with the solar 
magnetic activity. 



5.1. X-ray vs. stellar luminosities 

First, to know more about the X-ray properties of the 
members of the whole TTS sample (with upper limits if 
they are not detected with the HRI) , we examine whether 
a correlation exists between the X-ray luminosity and the 

^ We took 19.2' for the HRI field radius, because one of our 
sources (ROXRF31=SR9), is detected up to this angle from 
the axis in the Core F field. 

We have excluded three new ISOCAM Class II sources for 
which we have only the K magnitude, and thus no Av estimate. 

for these sources must be small, and/or Ay high, which 
implies a high upper limit on the intrinsic Lx. This does not 
affect the statistical results. 
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stellar luminosity (both corrected from extinction) , and if 
so, whether it is the same as the one found in p Oph by 
CMFA. We chose for each Core F X-ray source its lowest 
X-ray luminosity (including HRI upper limits) to mini- 
mize the effects of X-ray variability. For TTS undetected 
by the HRI, we estimate count rate upper limits (3.25 g)P^ 
To establish the existence of a linear correlation be- 
tween logLx and logL*, we performed three statistical 
tests using ASURV: Cox's proportional hazard model, the 
generalized Kendall r test, and Spearman's p test. The 
probability of the null hypothesis (i.e. that this corre- 
lation is not present) is < 10^* for each of the three 
tests. Thus, a strong linear correlation between logLx and 
logL* is indeed present. We found the linear regression 
coefficients by using the Estimation Maximization (EM) 
algorithm under Gaussian assumptions and the Buckley- 
James method. The Buckley- James method gave results 
similar to those of the EM algorithm, but with a larger 
uncertainty on the slope. As the Buckley- James method is 
semi-nonpar ametric, this suggests that the residuals of the 
linear correlation may be non-Gaussian. We thus conser- 
vatively keep the slope uncertainty given by the Buckley- 
James method. The Lx-L* correlation is then given by 
(seeFig.|): log(Lx/ergs-i) = (1.0±0.2) x \og{LjLQ) + 
30.1. We note that the censoring fraction is so high that 
the correlation line misses most of the data points and de- 
pends entirely on the location of the few lowest detections. 
The correlation dispersion may be due to X-ray variabil- 
ity, and also to TTS spectral type and age differences: 
Neuhauser et al. (1995) points out that the ratio Lx/i* 
increases with decreasing effective temperature, and shows 
a variation of Lx with age. 

This correlation spans three orders of magnitude in Lx 
and two in L^,. The slope of this correlation, a, is equal to 
1.0, and the TTS X-ray luminosity is then approximatively 
given by the simple proportionality: Lx/-^* ^ 10^"*. We 
thus confirm that the characteristic for TTS in p Oph 
is 10"''', with a large dispersion up to a level ^ lO^'^. 
There is no evidence for the "saturation" effect seen at this 
level in late type main sequence stars, and attributed to 
the complete filling of the stellar surface by active regions 
(Fleming et al. 1989). 

A similar correlation between Lx and L^, was found for 
Class II and Class III sources in the previous p Oph study 
of CMFA (the method used to estimate L^, was different, 
but very similar numerically; see j |3.2| ), but also in other 
star-forming regions: Chamaeleon (Feigelson et al. 1993; 
Lawson et al. 1996), Taurus- Auriga (Neuhauser et al. 

1996| ). However, the 
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1995), and IC 348 (Preibisch et al. 



slopes may not be identical: while Feigelson et al. ( 1993D 



When the YSO is in both Core A and Core F field, 
we use the longer Core F field exposure to estimate the 
count rate upper limit, using the EXSAS command COM- 
PUTE/UPPER_LIMIT S, and we use the extinction estimate 
from Bontemps et al. (200C) to compute the corresponding 
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Fig. 5. Intrinsic X-ray and stellar luminosity correlation 
for the "T Tauri star sample" . The squares represent the 
T Tauri stars detected with the HRI in Core A obser- 
vation, and the T Tauri stars detected with the HRI in 
Core F observations for which the lowest X-ray luminos- 
ity of the 3 observations is not an upper limit. Downward 
triangles correspond to the T Tauri stars detected with 
the HRI in Core F observations for which the lowest X- 
ray luminosity of the 3 observations is an upper limit. 
Arrows are the hundred upper limits for the HRI unde- 
tected T Tauri stars. The highest upper limits correspond 
to high extinction sources (see Fig. |^). The solid line shows 
the correlation found using ASURV: log(I/x/ergs^^) = 
(1.0±0.2) X log(i*/io) -I- 30.1. 



Preibisch et al. ( |1996D , and CMFA find the same a = 1 
slope as above. On the other hand, Lawson et al. (1996) 



limit on the intrinsic X-ray luminosity. 



found a = 0.55, on a better characterized, enlarged X-ray 
source sample in Chamaeleon, stessing the importance of 
having a sample as complete as possible. Nevertheless, the 
fact that we find the same slope as CMFA with an enlarged 
sample of X-ray sources in the same cloud is certainly a 
good internal consistency check between the PSPC and 
the HRI 

We find that the majority of the X-ray luminosity 
upper limits are above the Lx-Li, correlation. Only 5% 
of the X-ray undetected TTS are below the correlation, 
mixed with X-ray detected TTS. This is consistent with 
the idea that all TTS in p Oph may be X-ray emitters 
with ix/i* ^ lO"'^. Therefore, the TTS undetected by 
the HRI do not make up a separate population, but must 



10 



Grosso et al.: A combined ROSAT-HRI/near-to-mid IR study of the p Oph dark cloud 



have X-ray properties comparable to that of the detected 
population, verifying the same correlation. 

5.2. The X-ray undetected T Tauri star population 

Using the previous correlation between the stellar and X- 
ray luminosities, we can now estimate for each member 
of the TTS sample the intrinsic X-ray luminosity in the 
ROSAT energy band,0 and compare it with the HRI de- 
tection threshold to understand the X-ray detectability of 
the TTS sample with the HRI. However, the comparison 
is not straightforward, since the HRI detection threshold 
depends on both instrumental effects and extinction along 
the line of sight. 

Fig. ^ shows the instrumental effects: the HRI count 
rate threshold (3.25 a) increases away from the pointing 
axis. We interpret this dependence as the consequence 
of the point spread function degradation and reduced 
sensi tivity off-axis of the ROSAT mirrors (David et al. 
1997 ). With the X-ray spectrum assumptions described 



in §4.1, we have determined using EXSAS the conversion 
factor, /, between the HRI counts and the apparent X- 
ray luminosity (i.e., in the absence of extinction) in the 
ROSAT energy band (0.1"2.4keV), Lx.app- We find: / = 
6.8 X lO^^ergcts"^ s for d = 140 pc. The minimum X-ray 

ranges from 



7 X 10^^ ergs"^ 



luminosity for a 3.25 a HRI detection, Lx 
^ 7 X 10^''ergs~^ on-axis (angle=0') to 
off-axis (angle=19.2') (see Fig. ||). 

In case the X-ray sources suffer some extinction equiv- 
alent to Ay magnitudes, the values of ix.min on-axis and 
off-axis must be corrected to obtain the corresponding in- 
trinsic minimum X-ray luminosities as a function of Ay: 
if a source is heavily extinctcd, this minimum may be up 
to two orders of magnitude higher or more than in the 
absence of extinction (see for example the high values of 
the upper limits of Fig. |^). 

Fig. 1^ plots the X-ray luminosities of the TTS sam- 
ple as a function of (or A^h)- These X-ray luminosities 
were estimated from the stellar luminosities using the cor- 
relation discussed in the previous section. The points are 
compared with the ffi?/ threshold curves Lx.min = /(^v), 
computed both on-axis and off-axis. The HRI detected 
TTS (crossed dots) are found to be rather bright {^^ 
1029_ lO^^ergs"^) and weakly extincted (^y < 10). 
The undetected TTS have estimated X-ray luminosities 
below the computed HRI detection threshold. In particu- 
lar, we understand why the new ISOCAM Class II sources 



(Bontemps et al. 2000), characterized both by low stellar 
luminosities (^ 0.05 L0) (and thus presumably low pre- 
dicted X-ray luminosities, 6 x lO^^ergs"^), and rela- 



We note that this method attributes a larger X-ray lumi- 
nosity to the TTS having X-ray luminosities below the Lx-i* 
correlation. By this effect X-ray undetected TTS can be put 
above the HRI detection threshold, but this concerns only 5 
cases, see Fig. ^. 
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Fig. 6. HRI count rate thresholds (3.25 cr) vs. off-axis dis- 
tance {bottom scale). Left (resp. right) scale: HRI count 
rate (resp. apparent X-ray luminosity) threshold (3.25 cr). 
Plus signs (asterisks) represent Core A (F) T Tauri stars 
of the "T Tauri star sample" undetected with the HRI. 



lively high extinctions {Ay ~ 20), could not have been 
detected with our HRI observation.]^ 

Now ISOCAM cannot per se recognize Class III YSO 
among its sources without IR excess, but X-rays can. 
However, a reliable census of Class III sources in p Oph 
is de facto limited by the sensitivity of X-ray observa- 
tions: Fig. ^ shows that the number of detected Class III 
sources decreases for low Lx and high Ay; roughly speak- 
ing Class III sources are mainly detected above Lx,min ^ 
2^029.6 gj.gg-1 ^Qj. equivalently L^,_i„in=O.35L0), and below 

Ay, mill ^ 30. 

This strongly suggests that unknown Class III sources 
may exist. We can estimate their number at least in 
regions at the periphery of cloud cores, by using the 
fact that the WTTS/CTTS ratio (or equivalently the 
Class Ill/Class II source ratio) is ~ 1, and also that 
the HRI is equally sensitive to Class HI and Class II 
sources (see §4). In the ffi?///S'0Cj4M overlapping area, 
this ratio is 19/22 ~ 1; on a comparable area Martin et 
al. ( |1998| ) also found a WTTS/CTTS ratio ~ 1. Since 
the "TTS sample" comprises 88 Class II sources and 35 
Class HI sources above L^, ~ 0.03 Lq, we predict that 



^■^ A few faint sources were however detected in spite of being 
below the nominal HRI detection threshold: they were proba- 
bly in an X-ray flaring state at the time of the observations. 
In particular ROXRF32 = GY238, far below the HRI detec- 
tion threshold, was detected only in the third Core F exposure, 
which supports this interpretation. 
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Fig. 7. Detectability of the "T Tauri star sample" with the HRI vs. extinction. Right scale: stellar luminosity de- 
termined from near-IR photometry; left scale: X-ray intrinsic luminosity predicted from the Lx-Li, correlation (see 
Fig. H); top scale: extinction expressed in Nh; bottom scale: extinction expressed in Ay. "New" ISOCAM Class II 
sources are IR sources with IR excess discovered by ISOCAM. "New" Class III sources are IR sources characterized 
as YSO candidates from X-ray or radio observations, and for which /5'0 CAM observed no IR excess. Crosses indicate 
T Tauri stars detected with the HRI. Solid curves show the detection threshold for the ROSAT HRI (on- and off-axis), 
assuming a Raymond-Smith spectrum with fcTx^lkeV and an exposure time of 77182 s. Several T Tauri stars were 
detected despite being below the nominal HRI detection threshold: these sources have likely been detected in a flaring 
state (see text for details). Dashed-dotted curves show the XMM-EPIC (medium filter) detection threshold for the 
same spectrum and exposure time; we assumed a noise of 4 x 10^^ cts s^^ for the EPIC-pn, and deduced the detection 
threshold for the all EPIC instrument applying a 0.75 factor (see the XMM Users' Handbook; Dahlem et al. 1999). 
The dotted curves show the XMM-EPIC (medium filter) detection threshold for an high plasma temperature of 5 keV 
(flare). The X-ray detection of the "T Tauri star sample" is limited by the HRI sensitivity. XMM-EPIC thanks to 
its increased sensitivity and its energy range (0.2-12 keV), is less sensitive to extinction, and should be able to detect 
numerous new Class III sources, as well as the low-luminosity Class II sources discovered by ISOCAM 
(see text for details). 
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(88 X 19/22) - 35 ~ 40 Class III sources remain to be 
discovered in X-rays in the HRI/ISOCAM overlapping 
area above Lx 3 x 10^^ erg s^^. These sources are not 
seen now either (i) because they are too faint in X-rays 
[Lx < ix,min) "Or cquivalcntly from the existence of an 
Lx vs. Li, proportionality, too faint in stellar luminosity 
(L* < L*,mm)- or [a) too absorbed [Ay > ^v,min), or a 
combination of both. XMM-Newton will be an ideal tool 
to reveal such a large number of unknown Class III sources 
in the future (see iQ, but, as shown in the next section, 
we can already figure out their nature to a large extent. 

6. The unknown Class III source population 

In this section, we seek to characterize the suspected un- 
known Class III source population, which is likely to exist 
on the basis of the Class II source findings by ISOCAM 
and other IR observations. We will use (i) the spatial dis- 
tribution of all known Class II sources, and (ii) the extinc- 
tion map derived from C^^O observations. 

6.L Compared spatial distributions of the Class II and 
Class III sources 

The conventional wisdom is that Class III sources are de- 
scendants of Class II sources after dispersion of their disks 
(e.g., Lada 1987; AM). As cloud cores have an internal 
velocity dispersion, stars form with an initial mean veloc- 
ity distribution, implying that they drift away from their 
formation site (e.g., Feigelson 1996). This is the usual ex- 
planation for the increase with distance from cloud cores 
of the Class III/ Class II source ratio (or equivalently at 
that stage the WTTS/CTTS ratio), which is < 1 within 
the core region (e.g., CMFA), and reaches values 3> 1 far 
from the cores (e.g., Martin et al. 1998). This implies a 
larger spread of the spatial distribution of the Class III 
source population compared to that of the Class II source 
population. 

Let us first study the spatial distribution of the Class II 
source population within the HRI/ISOCAM area. We an- 
alyze the source surface density by using a 2-D Gaussian 
filter of a given FWHM on source position. The choice of 
the FWHM is optimized to enhance the contrast between 
regions of low and high source density, and thus reveal 
any clustering. Fig. ^ shows the resulting density map, in 
the form of dashed contours obtained with FWHM=6'. 
The Class II sources show three strong density peaks well 
centered on DCO"'' cores A, B and F, which is consistent 
with the idea that most of these sources were born in these 
cores. However, in spite of its comparable DCO''' line- 
of-sight density, core C appears much poorer in Class II 
sources; the weaker star-forming activity of this core is 
confirmed by th e pres ence of only one Class I source (see 
Bontemps et al. |2000| ). 

One can go one step further by comparing the source 
distribution with the matter distribution along the line- 



of-sight, i.e., with the extinction map. The 000+ radical 
is a good indicator of large densities (n ~ 10^-10^ cm~^) 
in cold cores, but the relevant regions occupy a relatively 
small volume; by contrast, C^^O, which is generally op- 
tically thin and sensitive to smaller densities, is a good 
column-density tracer. Using this molecule, Wilking & 
Lada ( 1983| ) derived an extinction map of the p Oph cloud 
center, showing that the denser regions have an equivalent 
visual extinction Ay between ~ 30 and ^100. 

Fig. ^ displays the C^^O contours, labeled in Ay by 
steps of Ay « 20, starting at Ay = 36, from Wilking 
& Lada ( 1983] ). Correspondingly, the Class II sources are 
represented by black dots of size decreasing with Ay , from 
low extinctions (Ay < 9) to high extinctions (Ay > 45), 
by steps of Ay « 10. A large majority of these sources 
are seen to have moderate extinctions (Ay < 18), even 
in the areas overlapping regions of high extinctions traced 
by C^*0. This implies that such sources are actually only 
moderately embedded in the cloud, in front of the densest 
regions traced by C^^O and DCO+, rather than within 
them. The spatial distribution of the Class II sources 
can thus be more appropriately described as gaussian-like 
three-dimensional overlapping "haloes" around the DCO"*" 
cores A, B and F. This also implies that at least a fraction 
of the Class II sources with high extinctions are not nec- 
essarily really embedded in the densest regions, but may 
be part of these haloes behind the dense cores. 

In the very same fashion. Fig. ^ displays the distri- 
bution of the known Class HI sources, this time using 
FWHM=8'j3 This distribution is different from that of 
Class II sources in Fig |[ With respect to the DCO"*" cores, 
there is a strong density peak ~ 5' SW of the location of 
core F, and no peak associated with cores A and B: in 
contrast with the distribution of Class II sources, there is 
a deficiency of Class HI sources in the cloud center regions 
with high extinction. 

The explanation for this apparent absence may be as 
follows. Whether they lie on the line-of-sight to regions 
of moderate extinction, or of high extinction, moderate- 
extinction Class II sources are found essentially every- 
where. Therefore we also expect to have low- extinction 
Class III sources everywhere, in a ^ 1:1 proportion. If they 
are not detected with the HRI, it can thus only mean that 
they are too faint in X-rays, hence have a small stellar 
luminosity. In addition, as for Class II sources, we must 
expect along the line-of-sight to the densest regions of the 
cloud to also have moderately embedded Class HI sources 
at the back of the cloud. Their spatial density should be 
roughly comparable to that of the unseen Class III sources 
in the front, i.e., yield a small absolute number given the 
compact size of the C^^O cores. 



A larger FWHM must be used because there are fewer 
Class III sources than Class II sources and because they are 
less clustered. 
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Fig. 8. Spatial distribution of Class II sources in the p 
Ophiuchi Cloud. Black dots show the position of Class II 
sources. The dot size gives information on the source vi- 
sual extinction. The dashed contour map is an estimate of 
the local Class II source surface density (in linear arbitray 
units) obtained using a sum of Gaussians (FWIIM=6') 
centered on each Class II source position; crosses show the 
maxima of these peaks. The letters show the location of 
DC0+(J=2-1) dense cores A, B, C, E, F (see Fig. |). The 
contour map shows the visual extinction (Av=36, 54, 72, 
90) d erived from C^^O column density (Wilking & Lada 
1983 ). A scale of 0.5 pc is shown for d = 140 pc. 




Fig. 9. Spatial distribution of Class III sources in the p 
Ophiuchi Cloud. White dots show the position of Class III 
sources excluding the early- type star SI. Crosses sh ow X- 



ray detected Class III source s (Montmerle e t al. 1983 



CMFA; Koyama et al. |1994 Kamata et al. |l997l ; this 
article). The dashed contour map is an estimate of the 
Class III source surface density (in arbitrary linear units) 
obtained using a sum of Gaussians (FWHM=8') centered 
on each Class III source position. Plus signs show the po- 
sitions of the density peaks of Class II sources (see Fig. ||). 
Class III sources with known spectral types are labeled, 
and put in an H-R diagram in Fig. QQ. 



There is also the possibility of having genuinely embed- 
ded (hence very young) Class III sources in these cores: we 
have no information about the Class Ill/Class II source 
ratio there, so it is impossible to estimate their number. 
Should this number be large (such that Class Ill/Class II 
> 1 for instance), this would be a problem for the earli- 
est stages of YSO evolution; one rather expects to have 
Class Ill/Class II ^ 1 if all stars form with a disk taking 
at least 10^ yr to dissipate. However, a disk stage is per- 
haps not necessary for very low-mass stars, which would 
increase the number of very young Class III sources. 

6. 2. Constraints on the nature of the unknown Class III 



Let us construct the H-R diagram of the 12 Class HI 
sources in HRI/ ISOCAM area for which we know the 
spectral types from the optical observations of Bouvier 
& Appenzeller ( |l992|) , and from the if-band observations 
of Luhman & Rieke ( 1999 ), using the stellar luminosities 



determined by Bontemps et al. ( 2000[ ).P^ Fig. [To| displays 
the result, along with the birthline and pre-main sequence 
evolutionary tracks of Palla & Stabler ( |1999| ) . According 
to these evolutionary tracks, the ages of the 12 Class III 
sources are found to be spread between ^ 0.2 and ^ 5 Myr. 

It is reasonable to assume that the unknown Class III 
sources have the same age spread. These sources are not 
yet detected in X-rays either because their intrinsic X- 
ray luminosities are too low, and/or because they have 
high extinctions. In the first case they have stellar lu- 
minosities below 0.35 Lq, and the isochrones imply that 
~ 0.6 M© for the oldest ones, and < 0.1 Mq for 
the youngest ones. In the second case, however, the un- 
known embedded Class HI sources can have luminosities 
higher than i*,min, i-e., so that they are not necessarily 
very low-mass stars. 

Unless the number of Class HI sources embedded in 
the densest regions is very high, our conclusion is that the 

The relative differences between these luminosities, and the 
bolometric corrected ones from spectral types are lower than 
20%. 
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Fig. 10. H-R diagram for the Class III source with a 
known spectral types. Dotted lines show Palla & Stahler 
( |1999D pre-main sequence tracks, continuous lines show 
isochrones. The bold continuous line represents the birth- 
line for the star build-up accretion rate (10~^ yr~^). 
White dots show the position of Class III sources with 
known spectral type (Bouvier & Appenzeller 1992; Luh- 
man & Rieke 1999; see also Greene & Meyer [l995 ) , and us- 
ing the Bontemps et al. ( ^000| ) luminosities. Crosses show 
X-ray detected Class III sources (Montmerle et al. 1983; 
CMFA; Koyama et al. |1994| ; Kamata et al. |1997| ; this arti- 
cle). We take d = 140 pc. The error bar shows the system- 
atic luminosity error corresponding to a distance error of 
±20 pc. 



hulk of the Class III sources which are undetected by the 
HRI and unrecognized by ISOC'AM should be made of very 
low-mass stars. 



7. Summary and conclusions 
7.1. Main observational results 

We have obtained two deep exposures of the p Oph cloud 
core region (d=140pc) with the ROSAT High Resolution 
Imager (core A: 51 ks, core F: 77 ks, in three partial ex- 
posures). The improved position accuracy (l"-6") with 
respect to previous recent X-ray observations {ROSAT 
PSPC, Casanova et al. 1995; and ASCA, Koyama et al. 



1994 and Kamata et al. 1997) have allowed us to re 



move a number of positional ambiguities for the detected 
sources. We have cross-correlated the X-ray positions with 
IR sources found in the ISOCAM survey of the same re- 
gion at 6.7 and 14.3 /im, in addition to sources known in 
the optical, IR, and radio from ground-based observations. 



We thus have now at our disposal the best-studied sample 
of X-ray emitting YSO in a star-forming region. We first 
summarize the main observational results of this article. 

(1) We detect 63 ffi?/ X-ray sources, and 55 are identified. 
Of the 55 identified X-ray sources 40 are PSPC sources, 
and 9 are ASCA sources. 

(2) The IR classification (ground-based and ISOCAM 
survey) for the 55 identified X-ray sources yields: one 
Class I protostar (YLW15=IRS43); 23 Class II sources, 
including 4 new ISOCAM Class II sources; 21 Class III 
sources, including 13 new Class III sources; 8 new 
Class II or Class III source candidates; one early- type 
Class III source (the young magnetic B3 star SI), and 
one field star (the F2V star HD148352). The contami- 
nation of the sample of new X-ray sources by field stars 
is negligible. 

(3) There is no statistically significant difference be- 
tween the X-ray luminosity functions of i7i?/-detected 
Class II and Class III sources, i.e. T Tauri stars with 
and without disks, confirming that the contribution of 
these disks to X-ray absorption, or emission (for in- 
stance by magnetic reconnection between the star and 
the disk), must be small. 

(4) X-ray variability of ffi?/-detected T Tauri stars has 
been studied by comparing the HRI data with the 
previously obtained PSPC data, and using HRI obser- 
vations done at three different epochs. The resulting 
statistics show that most of the sources are variable, 
and that their X-ray variability is consistent with a 
solar-like (hence magnetic) flare origin. 

(5) We use the information given both by the ISOCAM 
survey and by our HRI deep exposure to study the 
T Tauri star population of the p Oph dense cores. 
We confirm that essentially all Class II and Class III 
YSO are X-ray emitters, and that a strong correlation 
(log(Lx/ergs-i) = (1.0 ± 0.2) x log(L,/LQ) + 30.1) 
exists between the X-ray luminosity and the stellar lu- 
minosity of T Tauri stars, likely down to low luminosi- 
ties (L^ ^ 0.1 L0). We confirm that the characteristic 
Lx/i* for T Tauri stars is ~10~''' in the p Oph cloud, 
albeit with a large dispersion. There is no evidence 
for a magnetic "saturation" seen at a level of lO^'^ in 
late-type main sequence stars. 

(6) However, most of the new ISOCAM Class II sources 
are not detected by the HRI . We show that this is 
consistent with their intrinsic X-ray luminosities being 
too faint if "predicted" using the above Lx-L^, corre- 
lation. 



7.2. What have we learned ? 

(1) The first general conclusion we can draw from the HRI 
results presented above is a complete confirmation of 
the PSPC results obtained by CMFA. This was not a 
priori obvious, since the CMFA population {PSPC and 
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near-IR) overlaps, but is different from, the HRI/neai- 
IR/I^OC^Mpopulation presented in this paper: many 
PSPC sources are not detected by the HRI (see Ap- 
pendix C), and some HRI sources are Class II and 
Class III newly classified thanks to a combined iden- 
tification with ISOCAM. This shows that the Lx-Li, 
correlation is robust for the p Oph TTS. 
(2) The second, and perhaps most important, conclu- 
sion is the probable existence of ^ 40 unknown X- 
ray YSO down to a limit of Lx 3 x lO^^ergs"^ in 
the 7J7?///S'OCylM overlapping area, which should be 
mainly low- to very low-mass (< 0. 1-0.6 M©) diskless, 
"Class III TTS" . This prediction is based both on the 
use of the Lx^L^, correlation, legitimated by its robust- 
ness, and on the discovery of a large number of faint 
new IR sources by ISOCAM. As shown below, it may 
be soon verified by the next generation of X-ray satel- 
lites, namely XMM-Newton and Chandra. In this re- 
spect, the present paper can be taken as a "transition" 
paper between two generations of X-ray satellites. 

Why is the detection of these "unknown TTS" impor- 
tant ? Because they are diskless, they are unlikely to be 
recognized as YSO by IR observations alone; and because 
they are likely to be as numerous as the YSO with IR 
excess, they have to be included in any reliable census of 
YSO, with an impact on such basic quantities as the initial 
mass function, or the star formation efficiency, especially 
if considered from an evolutionary point of view. For in- 
stance, from the results in this paper it is impossible to 
study the real connection between the distributions of the 
Class II and Class III sources in the densest regions, in 
particular to see whether the distribution of the Class III 
sources is also centered on the same DCO"*" cores as the 
Class II sources. The number of Class III sources embed- 
ded in the densest regions may, or may not, be comparable 
to that of the Class II sources, depending on the timescale 
for disk dispersal, especially among low-mass YSO. An 
X-ray improved census of Class III sources may also be 
crucial in determining whether a burst of star formation 
is presently going on in the p Oph cores, as some recent 
indications suggest (see Martin et al. 1998). It will also 
allow to study the Lx-L* correlation for Class III and 
Class II sources seperately, which was not possible in this 
paper (§) due to insufficient statistics. 

7. 3. The potential of XMM-Newton and Chandra 

To quantify the prospects for improvement in the X-ray 
domain, we have computed the detection threshold for the 
X-ray camera EPIC aboard XMM-Newton, which was suc- 
cessfully launched in December 1999 (see Fig. |^). The im- 
proved sensitivity and enlarged energy range (0.5-12 keV) 
of EPIC will allow to detect the weak ISOCAM Class II 
sources, and also to discover numerous unknown faint or 
embedded Class III sources, in particular if they have high 



plasma temperatures (several keV) reached during flares, 
and extend the census of this population towards the low- 
mass end. In the best case, the XMM-Newton sensitivity 
will reach Lx ~ lO^^ergs""'^ for Ay ^ 20, for long expo- 
sures (>75ksec). This is nearly two orders of magnitude 
more sensitive than ROSAT. In case the faint Class III 
sources turn out to be so crowed that confusion problems 
arise, the excellent angular resolution of Chandra will be 
critical. 

In p Oph, there are already several identified bona 
fide a nd ca ndidate brown dwarfs (see review in Neuhauser 
et al. 1999 , and references therein), and four of them have 
been recently detected in X-ray using the ROSAT PSPC 
archive (Neuhauser et al. 



Neuhauser et al. have 
also shown that brown dwarfs could be X-ray emitters 
with the same ratio log(ix/i*) ^ —4 than for T Tauri 
stars. Thus Chandra and XMM-Newton should be able to 
detect many more of these objects with low stellar lumi- 
nosity and masses, shedding a new light on their nature 
and early evolution. 
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Appendix A: HRI X-ray source detection 

Source detec tion was done using EXSAS (Zimmer- 
mann et al. 1997 ), and the standard command DE- 
TECT/SOURCES. This command generates a local 
source detection by a sliding-window technique followed 
by a maximum likelihood test which compares the ob- 
served count distribution on the full resolution image 
(pixel of 0.5") to a model of the point spread function 
(PSF) and the local background (Cruddace et al. 1988 ). 
The "likelihood of existence" is defined as C = — IuPq: 
where Pq the probability of the null hypothesis that the 
observed distribution of counts is only due to a statisti- 
cal background fluctuation; L provides a maximum like- 
lihood measure for the presence of a source above the 
local background. We take £=6.8 as detection threshold 
(^0 < 0.0011; or > 3.25 cr for Gaussian statistics) as ar- 
gued in CMFA. 



For source detection, the HRI lepovt (David 1997) ad- 
vises to screen out lower and higher HRI Pulse Height An- 
alyzer (hereafter PHA) channels, which are found to have 
the highest background. However, source counts should al- 
ways be determined using all the 1-15 channels to cancel 
out the uneven HRI efficiency distribution across the de- 
tector area (S. Dobereiner, private communication). We 
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Fig. A.l. Sources designations from Table |A.l| - p3.l| . The image of Fig. |^ was smoothed using an edge detection 
algorithm. 



decided to search X-ray sources above a fixed detection 
threshold in channels 1-15 and 3-8.[^ EXSAS gave us list 

^® In the core F image #3, many 1-15 channel detections 
with high £ were not found in the 3-8 channel band Thus, 
we decided to take channels 2-8 instead of 3-8 in the three 
exposures. This third observation was the last of our program, 
one year after the Co re A o bservation (see Table |l|) . According 
to Prestwich et al. ( 199J ) , the mean pulse height decreases 
by 0.5 channels year" ; this effect might explains why we must 



of X-ray detections with positions, one sigma error box 
{<yx), ^, and count rate. We removed sources detected in 
channels 1-15 but not in 3-8, considering these detections 
as spurious. For instance, there are two hot spots in the 
Core A observation in the South-East corner of the HRI. 



decrease the lower channel boundary from 3 to 2. The upper 
channel boundary seems to be less sensitive, pr obably due to 
fewer counts in upper channels (see David et al. 1997). 
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Since hot spots are usually considered as spurious detec- 
tions, this criterion automatically removes them. 

The astrometry must be corrected from offsets of 
typically a few arcseconds due to the time-dependent 
boresight error in the ROSAT aspect system. To do this, 
IR or optical counterparts in a 10" radius circle around 
the X-ray sources were searched. We then selected a sam- 
ple of X-ray sources with an unambiguous counterpart 
and lax (half width) error box < 1", comparable to the 
IR/optical typical position error box (cr^ 1"). Then, 
offsets in right ascension {a) and declination {5) were 
estimated by individual offset weighted mean: Uoffset = 

SLl ^ ("X,i - a*,i), Soffset = Y17=l ^ ('^X,i - <5*a), 

with Wi = (l/cTxa + l/o-^.i)/ Z]r=i(l/'^x,i + l/o-^,i)- 
We found offsets ranging from -0.5" to 2.5". We 
subtracted these offsets, and checked the quality of 
our astrometry by estimating sample residuals mean 
before and after offsets subtraction: cr shift ~ n ^ 

I]"=l V(ax,i - - aoffsety + (<5x,i - - Soffset)'^- 

We found Ushift ranging from 0.5" to 1.2". Oshift and 
(7^ (= 1.2") were then quadratically added to ax to 
obtain an error box radius after astrometric correction 

i.'^totai = + ^.\»/t + <^l)- In the case of the three 
different core F observations, the images were aligned 
and merged after astrometric correction to obtain a 
single deep _ffiZ/ exposure of 77.2 ks. Source detection was 
subsequently performed as described in the article. 
Table 



A.l 



lists the HRI X-ray sources, for which we 
adopt, in Col. 1, the same acronym as in CMFA: "ROXR" 
(for p Oph X-ray ROSAT source), followed by "A" or 
"F" 1^ Fig. A.l indicates the source numbering. 

In order to allow easier comparisons with previous 
work. X-ray source positions are listed in both J2000 and 
B1950 equinoxes, with their latotai error box, in Cols. 2- 
6. The likelihood of existence C is in Col. 7. Count rates 
are indicated in Col. 8-11. For the core F field, the in- 
dicated positions [a, 6) and £ values correspond to ob- 
servation (#1, #2, or #3) where £ and the position ac- 
curacy are the best, i.e. when the count rate is highest. 
When an X-ray source is detected in one observation above 
the detection threshold, and not detected in other obser- 
vations, we have estimated the corresponding count rate 
upper hmits {3.25a), using the EXSAS command COM- 
PUTE/UPPER_LIMITS. We have noted that the detec- 
tion efficiency degrades with increasing angle to the axis, 
in the same way as the point spread function (this is dis- 
cussed in §5.2). 



The "legal" designation for a new ROSAT source is 
RXJHHMM.m ± DDMM, where J stands for coordinates in 
J2000 and m is the cut (not rounded) decimal value. Our ta- 
bles do not use this notation for simplicity in the discussion, 
but the correct designation is easy to reconstruct from the po- 
sition given if required. For example, R0XRF14 = ROXs20B 
= RXJ162714-2430. 



Appendix B: Optical/IR counterparts of the HRI 
X-ray sources 

We searched stellar counterparts for the 63 ROSAT HRI 
X-ray sources on the E SO/SERC second digitized sky sur- 
vey (DSS2). Fig. ^ gives the finding charts with BKLT 
IR sources for each of the 63 ROSAT Iffl/ X-ray sources. 



Table |B.l| gives identification lists for the two fields, 
and cross-identification with other surveys. Col. 1 is the 
ROXR numbering from detection (Table A.l). Cols. 2- 
4 are respectively cross-identification lists with the X-ray 



sources of CMFA (ROXRl), Casanova (|199|; R0XR2), 
and Kamata et al. ( |1997D . Dots mean "X-ray source un- 
detected", and dash "out of observation field" . Col. 5 gives 
the first name attributed to this counterpart. 

In the core A field, we find 26 X-ray sources, of 
which only one (ROXRAlO) remains without optical or IR 
counterpart. Of the 25 identified X-ray sources, 22 were 
seen with the ROSAT PSPC, and 4 are new detections 
(R0XRA3, 10, 16, 22). In the core F field, we find 37 X- 
ray sources, including 7 without optical or IR counterpart. 
Of the 30 identified X-ray sources, 18 were seen with the 
ROSAT PSPC, and 12 are new detections (R0XRF3, 8, 
12, 15, 18, 19, 24, 26, 28, 32, 35, 36). Altogether, 63 X- 
ray sources are detected, and 55 are identified. Of the 55 
identified X-ray sources 40 are PSPC sources. 

For sources with a low statistical significance (6.8 < 
C < 9.1, or IQ-^ <Po< 1-1 X 10-3; 3.25-3.9f7 for Gaus- 
sian statistics) we find X-ray sources with and without 
optical or IR counterparts. The X-ray sources without 
counterparts are always weak sources and may be spu- 
rious detections (locally high background), and this may 
therefore also be the case for weak X-ray sources with 
counterparts in case of chance spatial coincidence. For in- 
stance in the Core A field (respectively Core F) there are 
875 (resp. 1173) BKLT sources; this sample is dominated 
by background sources without detectable X-ray emission. 
To estimate the number of chance coincidences, we have 
placed in each field 10^ random X-ray source positions, 
and searched for each whether there is a BKLT source in a 
circle of 10" radius: we have then an estimate of the prob- 
ability to find a BKLT counterpart by chance within 10" 
from a spurious X-ray detection. This probability is 0.044 
(resp. 0.049) for Core A (resp. Core F), or approximately 
1/20 for both; in other words one (resp. two) spurious 
source identification are expected for the Core A (resp. 
Core F) field. As we have for Core A (resp. Core F) two 
(resp. 15) X-ray sources with £ < 9.1 out of 26 (resp. out 
of 37) , this implies that one weak X-ray source in Core A 
(resp. 13 in Core B) is real, which is consistent with the 
number of identifications. We are therefore confident that 
the identifications of weak X-ray sources with stellar coun- 
terparts are correct. 
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Table A.l. HRI X-ra.y sources in the p Oph cloud cores A & F. 
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Notes: ± gives the Idtotai error box; C is the hkehhood of existence, we give the maximum value for the observation 
set. 






Fig. B.l. Optical finding charts of the 63 ROSAT HRI X-ray sources of Tables |AJ and gj. Each map is a 30" x 30" 
extracted from the ESO/SERC sky survey red Schmidt plate using the second Digitized Sky Survey (one pixel=l"); 
North is at the top, East at the left. Circles show the ROSAT HRI 90% confidence error boxes (i.e. one sigma error 
box from Table A.l multiplied by 1.6). Asterisks show the BKLT infrared sources 90% confidence error boxes (~ 1.9"). 
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Fig. B.l. continued. 




Fig. B.l. continued. 
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Fig. B.l. continued. 



Appendix C: Comparison between HRI and PSPC 
observations 

Within the boundaries of our observation fields (core A 
and core F), there are 61 X-ray sources detected previ- 
ously with the ROSAT PSPC (53 from CMFA, and 8 
from Casanova 1994). However, 21 X-ray sources are not 
detected with the ROSAT HRI. This difference could be 
explained by lower observational sensitivity and/or source 
variability. To elucidate this point, we must estimate HRI 
count rates for PSPC sources, and compare them with the 
adopted 777?/ detection threshold (3.25 cr). 

We first estimate the conversion factor between the 
PSPC count rate in the energy range of CMFA (1.0- 
2.4 keV) and the HRI count rate in the whole energy range 
(0.1-2.4 keV). We did not select X-ray sources with am- 
biguous P SPC detection (10 sources with notes in Cols. 2- 
3 of Table B.l). For core F observation, the lowest detec- 
tion count rate was taken to minimize variability effects 
(four sources with upper limits are not selected): we kept 
only 26 X-ray sources. Since many of these sources are 
variable (as shown by the core F observations), a conver- 
sion factor estimator insensitive to extreme values of the 
sample is needed. This is why we take the median of the 
PSPC/ HRI count rate ratio, instead of the mean. We find 
PSPC i-2.4kcv count rate = 2.4 x HRI 0,1^2. 4kcV count 
rate. 



Fig. p.l| displays the HRI q i_2avcW count rate vs. the 
PSPC i-2.4kcv count rate. It shows two classes of sources: 
sources near the median, and sources beyond the median 
(with error bars) . The dispersion of points (within 1 rms) 
around the median value could be due to X-ray extinction 
effect on the conversion factor or to a variability factor 
< 2. Preibisch et al. (1996) calculated the conversion for 



the whole energy band of the ROSAT HRI assuming opti- 
cally thin plasma emission with Tx = 10^ K and different 
values for the X-ray extinction: for N^j increasing from 
6.5 X 10^^ cm~^ to 10^^ cm^^, the conversion factor de- 
creases from 2.5 to 2.0. Our observational estimate is in 
agreement with these values, which also show that the de- 
pendence of the conversion factor on X-ray extinction is 
small compared to the dispersion of count rates and can 
be neglected in our plot. We conclude that the dispersion 



is due to variability: WL20, GSS37, VSS27, and SR9 must 
have been in a high state during the PSPC observation, 
as were R0Xs4 and SR2 during the HRI observation, the 
other sources being essentially unchanged in both obser- 
vations. 

Using this conversion factor, we can estimate the HRI 
count rate from the PSPC count rate, and compare it with 
our HRI threshold computed with EXSAS. We find that 
14 sources are below this threshold (ROXRl-1, 6, 7, 16, 
20, 21, 27, 33, 34, 37, 40, 53, and ROXR2-16, 18), and 
the other 7 sources were in a high state during the PSPC 



slope=l/2.4 



o 




-1.0 



0.0 0.5 

log(PSPCi_ 



1.0 1.5 2.0 

^^,Y count ratc/ctsks"') 



Fig. C.l. Plot of ROSAT HRI (this paper) vs. PSPC 
(CMFA) count rates and sources variability. The dashed 
line is the median value of the conversion factor between 
HRI and PSPC count rates (2.4); dotted lines show the 
dispersion (Irms) around this median value. Count rate 
error boxes are shown for sources outward the dotted lines. 
Sources below (above) the dotted lines presumably flared 
during PSPC {HRI ) observations (see text for details). 
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ROXR logLx (d=140pc) 

ASCA NAME ISO IR Av 

1 2 Class [mag] [Lq] #1 #2 #3 #1+2+3 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

Al 2 — — SR2=SA0184375 — III* 0* 1.2* 30.4 

A2 3 — SR22 — II* 1.5* 0.4* 29.4 

A3 .. — B162538-242238 red nil 10.9 0.2 30.3 

A4 5* — IRS3 red nil 4.6 0.3 29.8 

A5 8 — 1 ROXs4=IRS10 blue nlll 5.4 1.7 31.4 

A6 9 — SR4 red II 1.9 1.5 29.9 

A7 10 — 2 GSS20 blue III 4.7 1.1 30.5 

A8 11 — Chini8*=SKS3 ... nlll* 0.5 0.2 29.1 

A9 12* — B162601-242945 blue nlll 8.1 0.5 30.2 

AlO .. — ? 

All 13 — 3 DoAr21=ROXs8 blue III 6.0 14.7 31.6 

A12 14 — VSSG19 blue nlll 3.9 0.5 29.9 

A13 15* — ROXCl=B162607-242742 blue nlll 20.6 0.7 31.3 

A14 18 — GSS29 red II 9.4 1.4 31.0 

A15 17 — DoAr24=ROXslOA red II 1.8 0.7 30.0 

A16 — . .* — WSB28 blue nlll 4.3 0.4 30.1 

A17 22* — GSS30-IRS2 blue III 19.5 0.9 30.9 

A18 24 — 4 GSS31 red II 6.1 5.9 30.4 

A19 25 — C3 S2 red II 11.8 3.7 30.4 

A20 26 — E124 red II 10.0 4.5 30.5 

A21 28 ~ 5 Sl=ROXsl4 bhTe m* 12^ 1100* 308 

A22 .. — B162636-241554 red nil 7.6 0.3 30.2 

A23 29 — C4 GSS37 red II 8.5 1.6 30.3 

A24 — 17 VSS27=ROXsl6 red II 5.6 1.9 30.7 

A25 32 — VSSG3 blue III 15.7 2.8 31.2 

A26 38 — SR21* red II 3.5 4.0 30.0 



Fl 23 — DoAr25* red II 0.7 0.8 29.7 <29.4 29.8 29.6 

F2 31 — GY112 blue nlll 3.6 0.4 <30.0 30.2 30.4 30.2 

F3 .. — GY122 blue nlll 2.7 0.08 <29.7 <29.6 29.8 29.6 

F4.. — < << ? 

F5 . . .... ? ? ? ? 

F6.. .... ? ? ? ? 

F7 35* — SR24S* red II 5.9 2.2 30.1 30.6 30.2 30.3 

F8 .. — GY193 blue nlll 7.4 0.3 <30.3 <30.2 30.0 30.0 

F9 36* — GY194 blue nlll 9.1 0.4 30.4 30.2 30.6 30.5 

FIO .. — < ^ <, l_ 

Fll . . — .... ? ? ? ? 

F12 . . — anonymous star* — II-III? < ^ ^ ? 

F13 .. — .... ? ? ? ? 

F14 — 22* — ROXs20B — III 2.7 0.3 30.3 30.4 30.4 30.4 

F15 .. — GY238 red nil 45.7 0.3 <32.9 <32.6 33.0 <32.5 

F16 39 — WL20 red II 16.5 0.9 <30.9 31.1 30.7 30.9 

F17 41 — 9A SR12AB=ROXs21 blue III 1.2 0.8 30.7 30.6 30.6 30.6 

F18 . . — B162720-243820 II-III? ? ? ? ? 

F19 .. — WSB49 ..* II* 1.9 0.2 <29.3 <29.1 29.3 <29.1 

F20 .. — ? ? ? ? 

F21 43* — YLW15=IRS43 red I 30* 10* 32.6 <31.1 <31.4 31.9 

F22 44 — VSSG25 red II 9.8 0.4 <30.8 31.7 30.4 31.3 

F23 46* — B162728-244803* — II-III? < < < ? 

F24 — . . — U0600-11613195* — II-III? ? ? ? ? 

F25 48* — BBRCG50 .... II-III? < < < ? 

F26 .. — B 162730-244726 blue nlll 10.3 0.6 30.6 30.5 <30.2 30.5 

F27 49 — GY292 red II 10.8 1.6 <30.6 30.7 30.8 30.7 

F28 .. — GY296* blue nlll 5.1 0.08 < < < 29.4 

F29 50 — IRS49 red II 10.1 1.1 31.6 30.4 30.4 30.9 

F30 51 — GY314 red II 6.4 0.8 30.4 30.4 30.3 30.4 

F31 52 ~ n SR9=ROXs29 red 11 IS 2978 29^2 29^9 29.7 

F32 . . — — GY377 blue nlll 16.0 0.4 < < < 30.6 

F33 54 — ROXs31=IRS55 blue III 6.0 1.8 30.7 30.7 30.7 30.7 

F34 55 — — WSB58 — II-III? ? ? ? ? 

F35 — . . — B162800-245340 — II-III? ? ? ? ? 

F36 — .. — B162812-245043 — II-III? ? ? ? ? 

F37 — 34 — HD148352 blue III* 29.9 30.0 29.8 29.9 



Comments to Table B.l 

•k = see notes below. ROXR A or ROXRF = X-ray source number (this articl e). R OXRl = CMFA X-ray source number. 
R0XR2 = Casanova (1994) X-ray source number. ASCA = Kamata et al. ( 1997 ) X-ray source number. Dash = out of 
observation field. Dots — unobserved source, red = ISOCAM source with IR excess, blue ~ ISOCAM source without 
IR excess, nil = new class II. nlll = new class III. ? — X-ray detected source for which intrinsic X-ray luminosity 
cannot be determined. < — X-ray undetected source. II-III? = class II or class III candidate (see Appendix B). 
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Notes to Table B.l 



R0XRA4 
R0XRA8 



ROXRAl: (7) The IR mdcx (aiR=dlog(AFA)/dlogA) is estnuated between 2. 18-4.69 from Walt er et ah ( |1994D and 
Jensen et ah (1997): we find am=-2.b. (8-9) These values are estimated from Walter et al. (1994). 
R0XRA2: (7) Martin et al. ( |1998D classify these stars as CTTS. (8-9) These values are estimated from GWAYL. 
(2) ROXRl-5 was identified with I RS3 a nd IRS5 (CMFA). 

(5) The position of Chini8 (Chini 1981) is 33" away from this optical star, the good position is given in 
Table 1 of Strom et al. ( |1995D by source number 3. (7) The IS0CAM\M2 and LW3 upper hmits exclude an IR excess 
(see Appendix D). 

R0XRA9: (2) ROXRl-12 was identified with an anonymous optical star (CMFA). 
R0XRA13: (2) ROXRl-15 was neither optical nor IR counterpart (CMFA). 

R0XRA16: (3) R OXR2 -16 was ~ 1' from R0XRA16 (identified with WSB28; CMFA), ie only two error boxes away, 
but Martin et al. (|l998|) identify ROXR2-16 with another star. 
R0XRA17: (2) ROXRl-22 was identified with GSS30-IRS1, 2, 3 (CM FA). 



1989|; A ndre et al. |l988| ). (8) The visual extinction 



R0XRA21: (7) SI is an embedded B-typ e star (Wilking et al. 
comes from the da ta of Lada & Wilking (1984), see Andre et al. ( 1988 ). (9) This B3-B5 stellar luminosity is taken 
from Andre et al. (|l988| ). 

ROXRA26: (5) In the BKLT survey this well known emission line star appears to be a 6.5" separation binary (this 
can also be suspected from the finding chart): the main component is B162710-241914 (J=8.56), the second one is 
B162710-241921 (J=11.27). 

ROXRFl: (5) The IR star B162623-244308 is also in the 90% confidence error box but its J-band luminosity (15.65) 
is lower than that of the well known emission line star DoAr25=B162623-244311 (9.29). 

R0XRF7: (2) ROXRl-22 was identified with SR24N and SR24S (CMFA). (5) In our observation #2 and #3 the 
counterpart of R0XRF7 is clearly associated with SR24S. In our observation #1 due to the weakness of the X-ray 
source the situation is less clear. We associate this source with SR24S. 
R0XRF9: (2) ROXRl-36 was ambiguously identified with GY193 and GY194 (CMFA). 

R0XRF12: (5) A weak star is visible in the 90% confidence error box on the DSS2 image, but this star is neither 
detected by the BKLT survey, nor by the PMM USNO-Al.O catalogue. 
R0XRF14: (3) ROXR2-22 was identified with ROXs20A and ROXs20B (CMFA). 

R0XRF19: (6) This source is just at the border of the ISOCAM survey. As a part of its flux is lost, Bontemps et al. 
(2000) do not characterize this source. (7) Martin et al. (1998) classify these stars as CTTS. We thus consider this 

source as class II. 

R0XRF21: (2) ROXRl-43 was identified with GY263 and IR S43 (CMFA). (8) Best value from Grosso et al. ( |l997|) . 
(9) This is the bolometric luminosity (see Grosso et al. 1997 ). 

ROXRF23: (2) ROXRl-46 was identified with an unnamed optical star (CMFA), probably the star B162730-244726. 
(5) We identified this X-ray with the IR star B162728-244803 just at the border o f the X-ray error box. 
ROXRF24: (5) This source is red in the PMM USNO-Al.O Catalogue (Monet et al. |l996D with B = 20.1 and R = 16.3. 
ROXRF25: (2) ROXRl-48 was identified with GY280, GY290, and GY291 (CMFA). 

ROXRF28: (5) The IR star B162735-244532:B is also in the 90% confidence error box but its J-band luminosity (> 17) 
is lower than GY296=B162735-244532:A (12.62). 

ROXRF37: (7) The Hipparcos distance is 75 pc: this star is a foreground F2V star. 



References to Table B.l 
( p^79| ). El 



B = Barsony et al. ( |1997D. BB RCG = Barsony et al. ( |1989|). D oAr = Dolidze & Arakelyan 
Elias (|1978| ). CSS = Grasdalen et al. (|1973| ). GY = Greene & Youn g (|1992] ). HD = The Henry Draper 



catalo gue (D raper |1918D . IRS = Wilking et al. ( |1989D . ROXC = Montm erle et al. (|T983^ROXs = Bouvier & Appen- 
zeller (^992|). S = abb reviat ion for "Source" in Grasdalen et al. ( |1973D . SKS = St rom e t al. ( |1995| Table 1)._SR_= 



Struve & Rudkj0bing(|l94|). U 
VSSG = Vrba et al. (|1975|). WL 



The PMM USNO-Al.O Catalogue (Monet et al. 199e). VSS 



Wilking & Lada (1983). WSB = Wilking et al. (1987). YLW 



= Vrba et al. (1976). 
Yoimg et al. (1986). 



observation (ROXRl-19, 30, 45, 47, and ROXR2-27, 30, 
33). We conclude that the non-detection of the 21 PSPC 
sources by the HRI can be fully explained by the difference 
in sensitivity and intrinsic variability. 

We can also compare our detections with ASCA 



(Koyama et al. 1994; Kamata et al. 1997; see observation 
field in Fig. l|) despite its lower angular resolution. Kamata 
et al. (1997) detected 19 X-ray sources, of which 10 were 



previously observed by Koyama et al. (1994). Compared 
to Einstein Observatory and i? T" FS'/'C observations, 
7 new X-ray sources were discovered by ^5(7^4]^ These 
7 X-ray sources are also not detected with the HRI. On 
the 12 sources already observed by Einstein Observatory 



Note that X-ray source 9B of Kamata et al. (1997) is in 
fact ROXRl-45 of CMFA. 
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and ROSAT PSPC, we detect 9 sources, the 3 others be- 
ing below our sensitivity threshold according to the above 
conversion factor. 



Appendix D: Optical/IR counterpart without IR 
classification 

Nine X-ray sources have optical/IR counterpart for which 
the IR classification is not known. Three of these X-ray 
sources are found in the ISOCAM survey, but with only 
upper limits in LW2 and LW3 filters. We g ive here their 
spectral energy distribution (see Fig. DJ), and discuss 
their possible IR classification. In case of doubt, the re- 
sulting Class II (or Class III) source candidates have not 
been included in the statistic studies of this article. 

ROXRA8: The counterpart of this X-ray source is the 
optical star Chini8=SKS3 (i? = 16.6, K = 9.5). The ISO- 
CAM LW2 and LW3 upper limits exclude an IR excess. 
We classify this source as a Class III source. 

ROXRri2: A weak star is visible in the 90% confi- 
dence error box on the DSS2 image, but this star is neither 
found in the BKLT survey, nor in the PMM USNO-Al.O 
catalogue. The low optical/near IR magnitudes imply a 
low luminosity for this object. We propose this source as 
a weak Class II or Class III source candidate detected dur- 
ing a strong X-ray flare. This source may also be a brown 
dwarf. 

ROXRF18: The counterpart of this X-ray source is 
the IR star B162720-243820 {K = 14.6). This star is vis- 
ible in the DSS2 (red) optical image, but it is not in the 
PMM USNO-Al.O catalogue (probably because only stars 
appearing both in blue and red images were accepted), 
thus we have no estimate of its B and R magnitudes. The 
low near IR magnitudes imply a low luminosity for this ob- 
ject. This source may be a weak Class II or Class III source 
candidate detected during a strong X-ray flare. ISOCAM 
LW2 and LW3 upper limits do not exclude an IR excess 
for this object. This object can also be a weak Class I 
protostar with a strong X-ray flare. 

ROXRF23: The counterpart of this X-ray source is 
the IR star B162728-244803 {K = 14.1; see note in Ta- 
ble B.l). The low near IR magnitudes imply a low lumi- 
nosity for this object. The SED of this source peaks in 
the i?-band. We propose this source as a weak Class II or 
Class III source candidate detected during a strong X-ray 
flare. 

ROXRF24: The counterpart of this X-ray source is 
an optical star {R = 16.3) named only in the PMM USNO- 



Al.O catalogue (Monet et al. 1996). Unfortunately, this ob- 



ject lies outside the BKLT survey. We propose this source 
as Class II or Class III source candidate. 

ROXRF25: The counterpart of this source is 
BBRCG50 observed only in K-hand {K = 11.7). This 
source is not retrieved in BKLT. The ISOCAM LW2 and 
LW3 upper limits exclude a strong IR excess. We propose 
this source as Class II or Class III source candidate. 
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Fig. D.l. Spectral Energy Distribution of optical/IR 
counterparts without IR classification. B and R data 
co me fro m the PMM USNO-Al.O catalogue (Monet et 
al. |1996|) ; J, H, K data come from BKLT; LW2 and 
LW3 are the ISOCAM upper hmit at lOmJy from Bon- 
temps et al. (200C). The dashed line shows the limit be- 
tween Class II source and Class III source classification 
(aiR=dlog(AFA)/dlogA=-1.5 from AM), "out" means that 
the source is outside the ISOCAM survey. 



ROXRF34: The counterpart of this X-ray source is 
the optical star WSB58=B162800-244819 {R = 17.3, K = 
9.3). Wilking et al. (1987) noted a probable Hq detection 
needing confirmation. The SED of this source peaks in the 
i?-band. We propose this source as Class II or Class III 
source candidate. 

ROXRF35: The counterpart of this X-ray source is 
the optical star B162800-245340 (R = 16.1, K = 9.6). The 
SED of this source peaks in the J-band. We propose this 
source as Class II or Class III source candidate. 

ROXRF36: The counterpart of this X-ray source is 
the optical star B162812-245043 (i? = 15.7, K = 9.4), 
which appears to be a close binar y ('-^ 1.5") in the sec- 



ond Digitized Sky Survey (see Fig. §^1). The SED of this 
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source peaks in the ff-band. We propose this source as 
Class II or Class III source candidate. 
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